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Abstract 

We investigate dynamical behavior of the equation of state of dark energy w& e by employing the 
linear-spline method in the region of low redshifts from observational data (Snla, BAO, CMB and 
12 H(z) data). The redshift is binned and Wd e is approximated by a linear expansion of redshift 
in each bin. We leave the divided points of redshift bins as free parameters of the model, the 
best-fitted values of divided points will represent the turning positions of Wd e where Wd e changes 
its evolving direction significantly (if there exist such turnings in our considered region). These 
turning points are natural divided points of redshift bins, and Wde between two nearby divided 
points can be well approximated by a linear expansion of redshift. We find two turning points of 
Wde in z € (0,1.8) and one turning point in z € (0,0.9), and u>d e (z) could be oscillating around 
w = —1. Moreover, we find that there is a 2a deviation of Wd e from —1 around z = 0.9 in both 
correlated and uncorrelated estimates. 
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I. INTRODUCTION 



It has been more than ten years since our universe was found to be in accelerating 
expansion A dominated and uniformly distributed energy component of the universe, 
called dark energyfDE), should be responsible for the acceleration. Many DE models have 
been proposed [2j-l6|. The simplest cosmological model is ACDM model, which contains a 
cosmological constant as dark energy. While ACDM model is still consistent well with all 
observational data, a lot of efforts have been made to find out whether DE is time-evolving 
or is just the cosmological constant. To do that, several parameterizations of equation of 
state (EoS) of DE have been proposed to fit with observational data, such as the ansatz 

n f] n 

Wde = + w ' z p]|j the EoS expanded by redshift, and the CPL parametrization [8|, |9[ 
Wde = wo + w a z/ (1 + z), expanded by scale factor. Both of them contain two free parameters: 
w , the present value of EoS, and w' or w a , represents the time evolution of EoS. Clearly, 
constraints of EoS obtained by using these parameterizations are model-dependent. Given 
an unreal assumption of EoS of DE, one ma y le ad to wrong conclusions. Some model- 



independent methods have also been proposed |10Ml3j] . such as the widely- used uncorrelated 
bandpower estimates (UBE) in which the redshift is binned and u^ e is assumed as a 

constant in each redshift bin. Note that the UBE method just approximates the actual Wd e 
by an averaged constant in each bin if DE is dynamical. If there are sufficient data, Wd e {z) 
can be accurately reconstructed. However, current data could only support a few bins, thus 
UBE is always used to test the deviation of Wd e from the cosmological constant and used as a 
supplementary for the parameterizations of Wd e - Note that the cubic-spline interpolation has 



also been proposed to study the binned Wde(z 
dynamic DE has been found 
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1g| . However, no convincing evidence of 
15l-ll7|. In addition, let us note that the ansatz, Wd e = wq + w'z, 



of redshift expansion and CPL parametrization exclude the possibility of an oscillation EoS, 
if they are used to fit the whole expansion history of the universe. While the UBE method 
needs enough bins to reveal the real dynamical behavior of DE, the errors will get larger as 
the number of bins increases. 

In this paper, we would like to probe the dynamical behavior of Wd e by using the linear- 
spline method. We will approximate Wd e in each redshift bin by a linear function w = 
Wq + w'z, and require that Wd e {z) is continuous in the region under consideration. Since 
most of data we used (e.g., Snla data) are in low redshift, we will focus on the region of low 
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redshift, such as z £ (0, 0.9) and z G (0, 1.8). In such regions, the width of each redshift bin 
is small and the linear expansion could be a better approximation of Wd e {z) than a constant 
in each bin. When fitting with the observational data, we leave the divided positions of bins 
Zi as free parameters. Since the linear function is monotonic, the best-fitted Z{ can represent 
the turning points of Wd e (z), where Wd e (z) is not linear enough or even non- monotonic (i.e., 
where d 2 Wd e /dz 2 departs from zero substantially). Actually we do find some turning points 
of Wde from observational data, and the constructed Wd e (z) just turns its evolution direction 
at the best-fitted positions of Zi. In this way, we only need to divide redshift into a few bins, 
the turning points are natural divided points of redshift and Wd e between two nearby points 
can be accurately reconstructed by linear expansion. Compared to the cubic-spline method, 
the linear-spline (LS) method can find turning locations of Wd e more accurately and reduce 
the errors due to less the number of bins. The LS method is also nearly model-independent, 
like the piecewise constant and cubic-spline method. Replacing the linear expansion by CPL 
parametrization in each bin, we have reached the almost same results. 

For the current status of observational data, LS method may be more suitable to study 
Wde than the piecewise constant and the cubic spline method. If DE is dynamical or even 
oscillating, by using the LS method it should be more possible to find deviations from the 
cosmological constant, at the turning points the deviation from —1 should be more explicit. 
If DE is just the cosmological constant, it seems more confident if the best-fitted linear 
expansions construct an w = — 1 line, while the oscillation of Wd e around w = —1 could 
disappear by averaging with the piecewise constant method. Compared to the piecewise 
constant case, the only price we pay is that there is one more parameter in the LS method if 
the number of bins is the same in two cases. Compared to the cubic spline method, the form 
of Wde(z) in each bin only depends on values of Wd e at two boundaries, thus the parameters 
in Wde{z) will not be heavily correlated. Furthermore, the cubic-spline method seems not 
suitable for finding the turning points of Wde- In all, the LS method could reconstruct Wde 
explicitly by using the least number of bins, and errors of the parameters from observational 
data will be small, compared to the case with more bins. 

The paper is organized as follows. In section II we introduce in detail the method we 
will use and construct corresponding cosmological models. In section III, we show how to fit 



18|], BAO data from SDSS DR7 [19j, CMB 



our model with 397 Constitution Snla sample 

datapoints (R, l a , z*) from WMAP5 2(3] and 12 Hubble evolution data 



21 



22|. The fitting 
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results are presented in section IV. We give our conclusions and discussions in section V. 



II. METHODOLOGY 

To fit models with observational data, we need to know the form of Hubble function H(z) 
(or E{z) = H(z)/H ). In a flat FRW universe 

e\z) = np)(i + zf + nf (i + z f + + zf + n^F(z), (i) 



where Q r , fl b , an< ^ ^efe are P resen t values of the dimensionless energy density for 
radiations, baryons, dark matter and dark energy, respectively, and Vti^ + + + 

)(0) 



f^ e = 1. The energy densities of baryons and dark matter are always written together as 
£lf\l + zf + + ^) 3 = fi^l + ^) 3 . The radiation density is the sum of photons and 



relativistic neutrinos 



20|: 



Qf) = Q(°)(l + 0.2271A^ n 
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where iV n is the number of neutrino species and = 2.469 x 10~ 5 /i~ 2 for T cm f, = 2.72hK 
(h = iifo/lOO Mpc • fcm • The evolving function F(z) for DE depends on Wd e {z): 

F(z) = e 3 f° 1 -^ d *. (2) 

For example, 

F(z) = (l + 2 ) 3 ( 1 +-o+- a ) e -^ ; 

for the CPL parametrization and F(z) = 1 for Wd e = — 1, respectively. Here we divide 
z G (0, 00) into m+1 bins and assume w^z) in the first m bins as 

W de (Zn~l < Z<Z n ) = W n ^ + w' n X (z - 2 n _ x ) , (1 < 71 < 77l) (3) 

and require Wd e (z) to be continuous at divided points: 

w n = uv-i + ^ x (z n - 2: n _i) , (1 < n < m - 1) (4) 

Note that here prime does not represent a derivative, instead is just the slope of the 
linear expansion in the n th bin. Thus the independent parameters are 

Wo,w' 1 ,w' 2 ,...,w' n ,...,w' m (5) 
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where the total number of parameters is 1 + m with m > 1. Alternatively, we can express 
w de as 

w ( z n) ~~ w ( z n-l) 

W de {z n _i < Z < Z n ) = w{Zn-i) H [Z-Zn-!), (1 < U < TTl) (6) 

z n Z n —i 

Now the parameters become w(z n )'s, which are values of Wde at the divided points and 
boundaries z n (0 < n < m) . In this case we have 

1 + z, 

x , ' ; 



where we have used zq = 0. For Wde in the last bin z G (z m , oo), we set it to be a constant 
wl, and 

F(z > z m ) = F{z m ){^^)^ +W ^ (8) 
i + z m 

Now the formula for H(z) is ready. 

There is one more thing to be mentioned: once we have fitted our model with the data 
introduced in the next section, errors of w(zi) are correlated, i.e., the errors of w(zi) are 
dependent on each other. New parameters can be defined by transforming the covariance 
matrix of w(zi), so that errors of new parameters are decorrelated and do not entangle with 



each other. The new uncorrelated parameters are referred to as the principal components 10, 



231 ]. and they are directly related to their own locations (unlike the correlated case). So errors 
of the uncorrelated parameters are more interpretable and meaningful. For more discussions 

nnn 

and implications of the uncorrelated parameters, we refer to the references pjj, H4|, 124J. In 
section IV, we will show both errors of correlated and uncorrelated parameters of Wde- The 
uncorrelated technique we adopt from [U] is as follows. 

1. Get the covariance matrix 

C = (WW T ) — (W}(W T ) (9) 

where W is the vector of w(zi). The Fisher matrix F is defined by F = C^ 1 . 

2. Diagonalize the Fisher matrix by an orthogonal matrix O 

F = T AO, (10) 
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where A is diagonal. 

3. Define a new matrix U as 

U = T A 1 / 2 0, (11) 

and normalize U so that the sum of its each row is equal to 1. 

4. Define new parameters g» by q = UW, where qi are components of the vector q. 
Clearly for the case of w(z) = —1 (i.e., the cosmological constant case), we will have 

qi = —1. The covariance of new parameters is 

(fa - fa»te - <*»> = ^^(^ < 12 > 

In this way the errors of the new parameters become uncorrelated. 

The uncorrelated parameters qi are linear combinations of w(zj), and the coefficients are 
just row elements of U. The transformation matrix U constructed in this method ensures 
that most of the coefficients are positive. So most of coefficients are in (0, 1). In this 
way, the original correlated parameters are weight-averaged, which leads to the uncorrelated 
parameters As a result, if Wd e is of the quintom form, the uncorrelated Wd e always looks 
more consistent with the cosmological constant than the correlated one. 



III. SETS OF OBSERVATIONAL DATA 



We will fit our model by employing some observational data including Snla, BAO, 
CMB and Hubble evolution data. The data for Snla are the 397 Constitution sam- 
ple f° r Snla is obtained by comparing theoretical distance modulus Hth( z ) = 
51og 10 [(l + z) J* dx/E{x)] + fio ( /iq = 42.384 — 51og 10 h ) with observed jiob of supernovae: 



X 



si I 



To reduce the effect of we expand xL with respect to Ho [251 ]: 

X 2 sn = A + 2Bno + Cnt (14) 



where 



[flthjZi] fig = 0) ~ Hob(Zi)] 

a 2 (zA 



2 



g _ Hthjzj] n — 0) - fi b(zj) _ y-> 1 



Eq. ()14p has a minimum as 



~2 _ 2 



A - B 2 /C 



which is independent of fiQ. In fact, it is equivalent to performing an uniform marginalization 
over /z , the difference between x 2 n an d the marginalized xt n 1S i us ^ a constant [25(. We will 
adopt x^n as the chi-square between theoretical model and Snla data. 



We will also use the Baryon Acoustic Oscillations (BAO) data from SDSS DR7 19], the 
datapoints are 

0.1390 ±0.0037 (16) 



r s {zd) 



and 



LV(0.275) 

£V(0.35) 
D v (0.2) 



1.736 ±0.065 



(17) 



where r s (zd) is the comoving sound horizon at the baryon drag epoch 26J], and 

1/3 



Dv(z) 



z 



z dx 
n H{x)J H{z) 



(18) 



encodes the visual distortion of a spherical object due to the non Euclidianity of a FRW 
space-time. 



The CMB datapoints we will use are (R,l a ,z^) from WMAP5 20] . z* is the redshift of 
recombination 28], R is the scaled distance to recombination 

dz 



R 



n 



(0) 



E(z) 

and l a is the angular scale of the sound horizon at recombination 

r(a#) 



(19) 



la = TT 



r s (a*) 



(20) 



where r(z) = J* dx/H(x) is the comoving distance and r s (a*) is the comoving sound horizon 
at recombination 

r a * n (n\ 1 

(21) 



r s {a^) = I - — - da, a t 



a 2 H(a) ' - 1 + z* 

where the sound speed c s (a) = 1/ \J 3(1 T Rbo) and Rb = 30^ /AVt^ is the photon-baryon 
energy density ratio. 

The x 2 °f th e CMB data is constructed as: 



2 _ vT r -l Y 
Xcmb — ^ U M ^ 



(22) 
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where 



X 



l a - 302.1 \ 

i? - 1.71 
z, - 1090.04 j 



(23) 



and the inverse covariance matrix 



n-l 



1.8 27.968 -1.103 
27.968 5667.577 -92.263 
-1.103 -92.263 2.923 



\ 



The fourth set of observational data is 12 Hubble evolution data from 



Xh is defined as 



12 



[H(zi) - H ob { Zi )f 



i=i 



at 



(24) 

2l| and |22j, its 

(25) 



Note that redshifts of these data fall in the region z G (0, 1.75). 
In summary, 



JCtotal "X.sn 4~ Xcmb Xbao ^CH' 



(26) 



IV. FITTING RESULTS 



A. Model I 

At first, we divide the whole region of redshift into 4 bins (i.e., m = 3), the divided 
points and boundaries are (0, z±, z%, 1.8, oo), where z\ and z 2 are left as free parameters of 
the model, and < z\ < z<± < 1.8. In the fourth bin we set wl = — 1. It means that we 
divide the region with z G (0, 1.8) into 3 bins and seek for two possible turning points of 
Wde(z) in this region. The reconstructed Wd e of the best-fitted model is shown in Fig. 



Models 


x 2 ■ 


X-cmb,min 


%-bao,min 


X,H,min 


%-total ,min 


Model I 


459.728 


0.204 


1.494 


5.983 


467.409 


CPL 


465.621 


0.251 


1.716 


10.818 


478.406 


ACDM 


465.759 


1.014 


1.470 


10.850 


479.093 



TABLE I: The best-fitted \ 2 of four data sets for Model I, CPL model and the ACDM model. 
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which indicates that there exist (at least) two turning points of Wd e in z G (0, 1.8) and the 
best-fitted values z x = 0.44 and z 2 = 1.07. Here Xtotai,min = 467.409 for the best-fitted 
parameters, roughly speaking it is a good improvement, compared with the corresponding 
X 2 t otai,min = 478.406 for the best-fitted CPL model and xlotai,min = 479.093 for the ACDM 
model. As shown in Table I, this improvement of Xtotai min * s mainly due to the decrease of 
Xsnmin anc ^ Xh min- ^ ^ s n °t surprising, as redshifts of the two data sets are distributed in 
whole range of z G (0, 1.8), while the BAO data are in the region z < 0.35 and the CMB 
data are in the region z G (0, z* ~ 1090). The result implies that these two data sets are 
quite favor of turnings of Wd e around z = 0.44 and z = 1.07 respectively. This result is 
consistent with recent UBE of Wd e |27J. While in CPL and ACDM models it is impossible 
to have such turnings of Wd e , which leads to the big differences between Xtotai min °f Model I 
and that of ACDM and CPL models. This result also implies that there exists the possibility 
with an oscillating EoS. Note that the error bar of Wd e in the third bin is larger than those 
in the first two bins because there are much less data points. 

We have also divided the region of z G (0, 1.8) into 4 bins, and found that there is almost 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



z 

FIG. 1: The best-fitted Wde for Model I (blue, solid), Model II (red, dashed) and the CPL model 
(green, dotted), the black line is for w = — 1. 
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no improvement of Xtotai min compared to the case of 3 bins, which indicates there is no more 
turning points of Wde in this region. 



TABLE II: The best-fitted parameters for Model I. 



parameters 



h zi z 2 u>(0) w(zi) w(z 2 ) u;(1.8) wl 



best-fitted values 



0.688 0.049 0.279 0.44 1.07 -0.63 -1.57 2.28 -16.84 {-1} 



B. Model II 

As data points with z > 1 are rather less than those with z < 1, to see clearly the 
evolution behavior of EoS in the region of low redshift, we now focus on the region with 
z G (0, 0.9), avoiding the possible turning point around z — 1, and set the divided points as: 
(0, z%, 0.9, oo), i.e., 



Wde\Z) 



w{0) + m(2i) 2 ; (0) z , 0<z< Zl 

w{zi) + "(O^frO fc _ Zl) t Zl < z < o.9 (27) 



-1 , 0.9 < z < oo 



In this case, we obtain the best-fitted tuning point Z\ = 0.45, and the best-fitted Wd e (z) is 
shown in Fig. [T] (the red, dashed line) which almost coincides with the best-fitted Wd e {z) 
of Model I in z G (0,0.9). This indicates that the data favor Wd e {z) to turn its evolution 
direction around z = 0.45, and favor an EoS with crossing the cosmoloeical constant (w = 

r~i 

— 1) [29]. Then we obtain la and 2a errors of parameters by using the MCMC method. 
Here we have fixed z\ = 0.45 in the process to obtain the errors of the parameters. Note 
that the errors of the parameters Wd e {zi) also represent errors of whole Wd e {z) in each bin, 
the la and 2a errors of Wd e {z) shown in Fig. 2 are obtained by connecting the corresponding 
error ranges of Wd e {zi). If another parameter set of Wd e (z) (as introduced in section II) was 
used, one will get the same result as that of Fig. 2. 

We see from the top left panel of Fig. [2] that there are deviations of Wd e from —1 around 
z = and z = 0.45 beyond la, and around z = 0.9 the deviation is beyond 2a. We 
decorrelate the parameters in Wd e {z) by using the technique introduced in section II. The 
uncorrelated errors are shown in the bottom left panel of Fig. [2j In that case, there are no 



10 



w L = - 1 




w L unfixed 




0.0 



0.2 



0.4 0.6 0.8 



-4 - 
0.0 



0.2 



0.4 



0.6 0.8 1.0 



0.5 



-0.5 



W L = - 1 

uncorrelated 




w L unfixed 
uncorrelated 




0.0 0.2 



0.4 0.6 0.8 



0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 2: lcr and 2a errors of Wde m Model II. Left panels are for the model with wl = —1 and right 
panels are with wl floating, top panels are for correlated parameters in Wd e { z ) an d the bottom 
panels are for uncorrelated ones. 

more explicit deviations from —1 around z = and z = 0.45, however, there is still a 2a 
deviation from —1 around z = 0.9. 

One may suspect that the explicit derivations are caused due to the fact that we have 
fixed the value of Wd e as wl = — 1 in the third bin. To check this, we consider the case with a 
floating wl. Two right panels of Fig. [2] show the correlated and uncorrelated results for the 
case with the floating wl- We see that in this case, there is even larger deviation from — 1 
around z = 0.9. We have also used the CPL parametrization to replace the linear expansion 
in each bin, and found that the errors are almost the same as those in the case of the linear 
expansion and there is still a deviation of Wd e from —1 around z = 0.9. We will extend 
our discussion of this result in the last section. In Fig. 3 we plot the likelihoods and weight 
functions of the uncorrelated parameters The weight functions are constructed from the 
transformation matrix U, which show how uncorrelated parameters qi are determined. It is 
shown that the parameters of Wd e are less correlated than that of the cubic-spline case 16]. 

We have also divided z G (0, 0.9) into three bins, to see whether there exist two turning 
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FIG. 3: Likelihoods and weight functions of uncorrelated parameters qi and corresponding weight 
functions. The left panels are for the model with wl = — 1 and the right panels are for the case 
with a floating wl- 



TABLE III: The best-fitted parameters and 2a errors for Model II with wl fixed to —1 or floating. 
ML is for "Maximum Likelihood", and the value in {} means this parameter has been fixed. 



parameters 


h 




ft (0) 


Z\ 


w(Q) 


w(zi) 


U>(0.9) 


W L 


best-fitted values 
ML and 2a errors 


0.684 

o.684±8:8g 


0.050 

50+ - 004 
u.uou_ .oq 3 


0.282 
0.286±g;8i 


0.45 
{0.45} 


-0.63 

u - D4 -0.38 


-1.59 
-1 fi4+ - 72 


1.24 

i 9 o+1.97 
i - zo -2.09 


{-1} 
{-1} 


best-fitted values 
ML and 2a errors 


0.687 

n fiS7+ - 027 
U.b8f_ _ 026 


0.049 
n n4q+ - 004 


0.280 

n noc+0.033 


{0.45} 
{0.45} 


-0.59 

u - ' -0.42 


-1.69 

+0.69 
L - ' u -0.79 


2.54 

q nq+2.58 
°- uo -3.66 


-1.72 
-1.771L 7 ? 



points of Wde in this region. We found that with the additional 2 parameters (z 2 and w(z 2 )), 
there is almost no improvement of Xmini compared to the 2 bins case (Model II). This 
indicates that there is no more turning points and Wd e {z) can be well approximated by just 
two linear expansions in the region z G (0,0.9). Of course, there is another possibility that 
the current data are not enough to find out more turning points. 
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V. CONCLUSIONS AND DISCUSSIONS 



In this paper we have investigated the dynamical behavior of the EoS of DE in the region 
of low redshift in a nearly model-independent way. The redshift in that region is binned and 
Wde in each bin is approximated by a linear expansion of redshift 2, and in the large redshift 
region we set Wd e to be a constant wl. While fitting the model with some observational data 
which include Snla, BAO, CMB and Hubble evolution data, we leave the divided points of 
bins as free parameters. If the evolution of Wd e is not monotonous, or is not linear enough 
in the region under consideration, the best-fitted divided points will represent the turning 
points, where Wd e changes its evolving direction significantly. In this way we can explicitly 
reconstruct Wd e by using a few bins, and the errors of parameters from observational data 
will be small due to the small number of bins. First we have tried to find the turning points 
within the region of redshift z G (0, 1.8), and set Wl = —1 in the region z G (1.8, 00) (Model 
I). Our results show that the data favor two turning points of Wd e in z G (0, 1.8), and Wd e 
may have an oscillation form 30(. Our results are consistent with those by the UBE method 



111 



27|. 



Since the main data points are in z G (0, 1) and our result in Model I shows there may be 
a turning point around z ~ 1, to see clearly the dynamical behavior of EoS in that region, 
we have focused on the region z G (0, 0.9) in Model II. We have found one turning point 
only in z G (0, 0.9), the reconstructed Wd e in the best-fitted model is almost the same as that 
reconstructed in Model I in z G (0, 0.9). We have also obtained the errors of Wde at la and 
2a in z G (0, 0.9). In both correlated and uncorrelated estimates with a fixed wl = — 1 or a 
floating constant Wl, we found that there is a 2a deviation of Wd e from —1 around z = 0.9. 

It is interesting to see whether the deviation of EoS from —1 around z = 0.9 is physical, 
or is caused by some unknown technical causes in fitting. If it is physical, it then clearly 
shows that DE is dynamical. But in UBE of Wd e there seems no such distinct deviation 
around z = 0.9, it may be due to the difference between the discontinuity of Wd e in the 



piecewise constant case and the continuity in LS case [3l|]. In [16], where the cubic-spline 
method is used, there is also no such an explicit deviation around z = 0.9, but it is likely due 
to its set of EoS in the last bin wl = w(l): to fit well with the data of z > 1, w(l) should 
be much minus, which would suppress the reconstructed Wd e around z = 1. Of course, it 
is also possible that such a big deviation around z = 0.9 is due to the non-smoothness of 
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Wde at the divided points in our LS method. In the LS method, Wde is continuous, but not 
smooth at the divided points, i.e., its derivative is not continuous at those points. In fact, 
Wde in LS method can be made smooth at the divided points, such as by the relation 

mil i / z— zi—-\ \ 

, \ sr^ w i - , \ , coshf — ^— -) 

= + E + Ah ooshtWAV ' <28) 



1=1 

■th 



where w\ is the slope of linear expansion in the i bin (i > 1 and W7 = 0), and A is related 
to the smoothed extent at the divided points. With this parameterization, one can still find 
out the turning positions of Wd e that are favored by observational data, and perturbations 
of DE can be calculated. 



Furthermore our results are also dependent on the data set we have used 32|. For 



example, although there is still a 2a deviation at z — 0.9 by using the widely-used data 



set Snla + CMB-shift R + BAO parameter A 33], now the best-fitted turning point in 
z G (0, 0.9) changes to z — 0.39. Whatever, from the observational data we have used, a big 
deviation of Wd e from —1 around z = 0.9 is found. Unlike the deviations around z = and 
z = 0.45, this deviation around z = 0.9 does not to be reduced in the uncorrelated estimate. 
At least, our result shows that if DE is dynamical it is more possible to find the deviation 
of Wde from —1 around this redshift value. 

If the EoS of DE is indeed of an oscillating behavior around —1, it is then not surprising 
that the cosmological constant always fits well with observational data because the oscillating 
behavior could be smeared in the luminosity distance. However, if the oscillation region of 
Wde is wide enough (like the case of our best-fitted Wde in Model I), DE may be distinguished 
confidently from the cosmological constant by more precise astronomical observations in the 
next generation. In addition, let us mention that an oscillating behavior of Wd e is also 
possibly due to some systematic errors in observational data, or due to some interactions 



between DE and dark matter 34 1. 
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